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S ■ Abstract 

In this paper we study an impurity in Bose-Einstein-Condensate sys- 
tem at T = OK and suppose the contact forms for boson-boson and 
boson-impurity interactions. Using Bogoliubov theory and a further ap- 
proximation corresponding to only think over the forward scattering of 
impurity by bosons, we derive a reduced Hamiltonian whose form is the 
same as the Frohlich Hamiltonian for large polaron. By using Lee-Low- 
qq \ Pines (LLP) theory for large polaron, we obtain the effective mass of 

impurity, the phonon number carried by impurity and the energy re- 
lated to the existence of impurity. In addition, we also discuss the valid 
range for forward-scattering approximation. 
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1 Introduction 

Since the impressive experimental achievement of Bose-Einstein-Condensation (BEC) [H0 
[3] quantum degenerate atomic gases have become one of the hottest domains. Temperature 
in such systems can be thought to approach zero so that the thermal fluctuations and 
classical phase transitions associated with thermal fluctuations are fully suppressed. But 
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quantum fluctuations due to Heisenberg uncertainty principle still exist and may be strong 
enough to lead to a quantum phase transition. So these systems are good playground to 
realization of quantum phase transitions [H El E] • Except probing for fundamental physics 
in the realm of dilute atomic quantum gas, BEC system is also an auxiliary but effective 
tool to sympathetically cool some systems, which can not be cooled by traditional cooling 
technique, to degenerate temperature [H [9]. So a quantum degenerate mixing system 
could be obtained. Extending such mixing system to a extreme case that there is only one 
atom belonging to sympathetically cooled system, we will expect a BEC system with an 
impurity. Fortunately, such system has been able to create in experiment by many ways 
[9| 1 101 \TT\ I12j. An impurity in BEC is very similar to polaron, an electron in ionic crystal. 
When the electron moves, it polarizes the lattice around it and inversely experiences an 
effective potential from the lattice. And in the BEC system, the existence of an impurity 
would modify the distribution of atoms in BEC and also produce an effective potential to 
act on itself [H] . In [T31 El CE] > their authors made good use of the product wave function 
similar to Landau-Pekar description to polaron |16[ [T7] to deal with the self-localization of 
impurity. The product wave function, the essence of which is the separation of freedoms, is 
a good approximation for polaron owing to large mass ratio of ion to the electron to make 
the Born-Oppenheimer approximation be valid. But the validity of product wave function 
is questionable in the BEC system with an impurity, owing to similar masses between the 
bosons and the impurity, and some deviations from it can be expected |15| . 

In this article, we start from the bosonic system Hamiltonian with an impurity, and 
by using Bogoliubov theory, derive an equivalent Hamiltonian which is the same as the 
Frohlich Hamiltonian [18119] for large polaron in section 2. In this equivalent Hamiltonian, 
Bogoliubov phonon plays a role of the phonon which comes from the lattice oscillation. Then 
we use intermediate coupling theory [20J of Lee-Low-Pines (LLP) for polaron to calculate 
the effective mass of impurity, the phonon number carried by impurity and the impurity 
energy in section 3. Finally, some discussions about validity of our approximate decision on 
condensate fraction and conclusions are given in section 4. 



2 Equivalent Hamiltonian for the Bosonic System with an 
Impurity 

We consider a single impurity immersed in a homogeneous BEC in three dimension space. 
Assuming that the boson-boson and boson-impurity interactions can be described by contact 
interactions, the many-body Hamiltonian [T3| reads 

where tub, mj are the masses of boson and impurity respectively, g = Ai:fi 2 aBB /fnBi 
A = 2-Kh 2 aiB /ni rec i are the interaction strength of contact potentials for boson-boson and 
boson-impurity interaction with the reduced mass m re d = mims / {rni + ms). fi, f rep- 
resent the coordinates of bosons and impurity, respectively. In the second quantization 
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representation, we can expand the bosonic field operator on a set of plane wave basis and 
change Hamiltonian ([1]) into 



k,p,q 

+A^4a ? exp -i{k-q)-f +— (2) 



where = h 2 k 2 /(2m,B) and P is the momentum of impurity. The chemical potential 
ix is introduced to keep the average bosonic number to be a constant. ak,a\. are the 
bosonic annihilation and creation operators conforming to the canonical commutation re- 
lation [ajfejOl] = Sf. q . Owing to only one impurity, the statistical property of impurity is 
unimportant and the same results are obtained for the fermionic and bosonic impurity. For 
the convenience, we suppose our system to be unit volume. 

At T = OK, the bosons macroscopically occupy the lowest energy state with wave vector 
k = 0, which makes it possible to neglect quantum property of ao, oj ana ^ consider them to 



be classical number [21]. Following the general procedure, we decompose at as follows 

a-fc ->■ \/no$ko + Ofc (3) 

no is the number density of atom on the lowest energy state with no = a/jOo- Substituting 
([3]) into ([2]) and only keeping terms until the second order, the Hamiltonian reads 

P 2 9 

H = - /xn + 2 n o + An + \fno{-^ + 5™o + A)(a + %) 

+A- N /np , y^ \ak exp(zfc • f) + a\exp(-ik • r) + a\a q exp[— - • f] 

+ H( £ fc - z 2 ) !^ + H( a ^ a -^ + 4 4°fe + a W-k) ( 4 ) 

where the sum with character J2k represents exclusion for k = state. Note that both 
interactions of boson-boson and boson-impurity make the BEC deplete. We take into con- 
sideration contributions of such two kinds of interaction to the condensate depletion at the 
same time by asking for the coefficient of ao,c4 t° vanish [21] and find 

li = gn + \ (5) 

When BEC happens, the gauge symmetry of the system is broken and a long-wavelength 
gapless mode exists which is imposed by Goldstone theorem |22j . In order to satisfy above 
condition, we find that if we make below approximation for boson-impurity scattering term 

A E 4 a « ex p [-^ - 9") • r i ~ A E a i afc ( 6 ) 

k,q k 

the gapless property of low energy excitations is met as we could see below. Physically 
this approximation implies that when the impurity is scattered by bosons, the momentum 
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of bosons and impurity will not alter. In other words, impurity only experiences forward 
scattering. So the Hamiltonian is further reduced to 



H = 



2mi 2 u 2 / 



+A- v /rap \a k exp {ik • f) + a\ exp {—ik • r) 
k 

+\ 5Z( £fc + 9no){a\a k + a^ k a ] _ k ) 
1 k 

+gn (a k a- k + a\a ] _ k ) (7) 
Making the Bogoliubov transformation 

v-kh - v k tf_ 



a 



ak = u k b k - v k b_ k 
t 



-v k b k + u k b ] _ k (8) 



and imposing the conditions 



We obtain 



2 2 1 

u k -v k = l 

i u k + v l)9no ~ 2u k v k (gn + e k ) = 

£fc = («fc + ' y fe)( £ ik + gno) ~ 2u k v k gn (9) 



F = E B + ^- + Y,^b{b k 

k 

t 

+ J2[ V kbk^p(ik-f) + V k *ble^(-ik-r)\ (10) 



where V k = X^/n^(u k - v k ) and E B = -§ + ± Efc(£fc ~ e fc - 3^0 )■ £fc = \f^r2gn^e k 
with u| = u| — 1 = l( gfc±p° — 1), UfcWfc = It is very clear that low energy excitation ^ 
is gapless and completely same as that in the situation without the impurity, no is decided 
by atomic number conservation 

n = n + ^2 < a\a k > 

k 

= no + $^Vfc + $^(ufc + vjfe) <&£&*> ~^u k v k < &J&L fc + &*&-* > (11) 

k k k 

Where < • • • > represents the ensemble average. In the pure bosonic system, Bogoliubov 
transformation makes the Hamiltonian of bosonic system diagonal about phonon operators, 
so the last two terms on the right hand side disappear and no is decided by equation 

n = n +J24 (12) 

k 
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While for the bosonic system with an impurity the last two terms are nonzero, leading to 
that the particle number equation and bosonic condensate fraction are modified in contrast 
to pure bosonic system. In fact, the last two terms in (|lip illustrate the contribution of 
impurity to condensate depletion. But as we will see below, the last two terms are much 
smaller than other terms in (jllj) and are negligible due to small ratio of impurity number to 
boson number. Thus, we determine condensate fraction according to (|12p . in other words, 
we do not consider the contribution of boson-impurity interaction to condensate depletion. 
When no is determined, the Hamiltonian (jlOp is completely decided. 

Owing to nonexistence of zero-momentum phonon, exclusion for k = state in (jlOf) is not 
important and we can neglect this exclusion. In addition to a constant term Eb, the form of 
Hamiltonian (110p is the same as that of Frdhlich large polaron theory. Bogoliubov phonon is 
akin to phonon coming from lattice oscillations. When an electron moves in ionic crystal, it 
will be influenced more strongly by optical phonon, which represents relative movement of 
positive and negative ions and is accompanied by polarized electric field, than acoustic one 
which represents the movement of mass center and can not produce polarized electric field. 
So in usual Frdhlich large polaron theory, the phonon is considered to be optical mode one 
and its frequency is chosen to be independent of the wavevector and be a nonzero constant 
due to finite energy gap for optical phonon. But in BEC system with an impurity phonon 
is gapless, moreover we must consider the dependence of dispersion relation on wavevector. 
As a summarization of this section, we get an effective Hamiltonian similar to Frdhlich large 
polaron by making Bogoliubov and forward-scattering approximations. 



3 The Effective Mass, Phonon Number and Impurity Energy 

Below, we follow the LLP approach [20] to calculate some properties of such system. The 
essence of the LLP theory consists in combining the canonical transformation with the 
variational principle to get approximate ground state of the system. Supposing |<1? > is the 
ground state H\& >= E\& >, we make the canonical transformation |3> >= U\\^f > with 
U\ = exp [—iJ2k b k b k k • r\ and the Hamiltonian is transformed into H = U{~ HU\ 

H > = Ehl (P- E fc hkbjb k f 

+ E ^b{h + J^(V k b k + V k *b{) (13) 

k k 

In (|13p . the coordinate of impurity is expunged and its momentum becomes a good quantum 
number. In fact, the momentum P in (|13p represents the whole momentum of the system. 
Introducing the variational function f(k) and making another canonical transformation 
>= C/2 1 > with U2 = exp \}2 k {b\f{k) — b k f*(k))] and |0 > representing the vacuum of 
the phonon, the Hamiltonian (|13|) can be rewritten as H = U^H Ui 

s = e b + {p - E 2 k *f kh)2 + E [Vkf(k) + v k *nk)] + ^ e k • q v(k)\ 2 \m\ 2 
+ E i/( fc )i 2 L-—P-k + d + E tkblh + E b l b k E — * • $f(<i)\ 2 



mi 



mi 
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+ E *: 



V£ + f(k) U k + e k P . k + —Y J k-q\f{q)\ 



mi 



mi 



+ h.c. 



+ 2m~ r E 1 ■ q\blblf(k)f(q) + 2b\b q f(k)r{q) + W (*)/* (?) 



k,q 



+ —J2k.q\f*(q)blb k b q + f(q)b\blb k 



mj 



k.q 
2 7„2 



(14) 



where e k = h k /(2m/). The ground state energy of the system E q (P) =< 0|-ff|0 >= 
£ B + £(P), where 



£(P) 



2m/ 



+ 



2m/ 



+ E^/( A; ) + ^/*( A; )] + Ei/( fe )i 5 



E^i/( fc )P 



mi 



(15) 



and f(k) is decided by = = 

V k + f*(k) 



Ck + e k -—P-k + —Y,\f(Q)\ 2 k-q 
mi mi q 



According to the LLP approach, letting J2k ^k \ f(k) | 2 = rjP, we have 



f*(k) 



^k + e k -^(l- V )P-k 



(16) 



(17) 



After introducing the parameter ij, the behavior of the system completely depends on the 
value of m We further have 



,p = E — 

k [& + €t-i(l-l/)P^ 



(18) 



The parameter rj is self-consistently decided in (|18p . which is the main equation in LLP 
theory. Substituting (|17|) into the expression of E(P) 



E(P) 



2mi 



a-^ 2 )-E 



t^k + e k -^(l-r,)P-k 



(19) 



Below we make a small quantity expansion about P to (fl~8l) till the second order to get 



(20) 
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Choosing the orientation of P to be along the z axis, we have 

12 o*2 



\ " 



\Vk\ 2 2ft' 2 2 

rrr K COS 



(21) 



In order to calculating the effective mass of impurity, we also expand the E(P) to second 
order of P and find 



Vk\' z | P 2 
^ 6 + £fc 2m/ 



W = -T.T^ + ^-V 2 ) (22) 



P2 (i _ v f V J^ 1 ' — A: 2 cos 2 9 



2m 7 ^ (£ fc + e k f mi 

By means of (I2ip . the effective mass of the impurity m 7 

m 7 = 23) 

1 — 7) 

Introducing another parameter W = ??/(l — ff) which, in fact, can reflect the effective 
strength of coupling in contrast to the situation for polaron, the effective mass can be 
expressed into mj/mi = 1 + W . Integrating the angular variable, W is reduced as 

8p m B A + m B \ 2 (k B ai B ) 2 n ( ) 



W = — — 1 + —2. v^±Bl n w (24) 

3V3n^ mj \ ^ mj J (k B a BB yl 2 1 ' 



where 



/•CO 

n(°) = / dy 
Jo 



x ^2 



(y 2 + 2p)V2[(y2 + 2p )l/2 + (mg/m^p 

(25) 



Where in order to scale the expression we define a parameter k B like the Fermi wavevector 
n = k B /(6TT 2 ). In addition, we also introduce another notation p = no/n for condensate 
fraction. From the expression of W and under the approximation (|12p which signifies the 
condensate fraction p is independent of aj B , we easily find that W is proportional to the 
square of the ai B for definite the a BB and m B /mi. Fig. 1(a) apparently respects this kind of 
behavior. The dependence of W on m B /mi for definite a BB and ai B appears to be indirect 
due to the dependence of D^ ) on m B /mi. Fig. 1(b) show the corresponding change of the 
effective coupling strength as function of m B /mi. 

At zero temperature, the number of phonon is zero for purely BEC while there exists 
some phonon excitations in BEC with an impurity. Below we calculate the phonon number 
carried by the impurity n ph = £ fc < >= Efc \f( k )\ 2 

2p (k B a IB ) 2 ( m B \ 2 (1) 

n p h = r- o /9 Ti u/2 1 + 11 v 26 ) 

V37T 3 / 2 (k B a BB yi 2 V mi J 

where 



nW = / dy - - — (27) 

Jo (y 2 + 2p) 1 /2 [( y 2 + 2p )i/2 + (m B / mi )yY - 4(m B /m 7 ) 2 (l - r]) 2 P 2 
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(a) 



3 |B k B 



(b) 



m B /m, 



Figure 1: The behavior of the effective mass m}/mj = 1 + W as a function of clib for 
different mass ratios tub /mi of atom to impurity (a); as a function of mass ratio mg/mj 
for different boson-impurity scattering length clib (b). We have chosen clbb^b = 0.87. 

and P = P / (2mBng) 1 ' 2 . When an impurity moves in BEC, the state of this impu- 
rity is sensitively dependent on its velocity according to superfluid Landau theory [23] . 
Landau used simple kinematic arguments to derive an expression for the critical velocity 
Vl = min[C(k)/k] where C(k) is the energy of an elementary excitation with momentum k. 
If the velocity of impurity is larger than Vl, the impurity will dissipate energy by colliding 
with BEC, but if not, the impurity will move without dissipation. Such behavior has been 
predicted theoretically [24] and observed experimentally |10j : as the impurity velocity is 
diminished below Landau critical velocity, there exists a dramatic reduction in the proba- 
bility of collisions. Recalling that in the process of deriving effective Hamiltonian (jlCJ we 
have made forward-scattering approximation. Only forward scattering corresponds to no 
scattering and vanishing collision probability. Thus we think that the reduction of collision 
probability for low impurity speed smaller than Landau critical velocity can guarantee the 
validity of forward-scattering approximation. In order to keep this approximation valid, we 
calculate phonon number n p h at P = 0. In ionic crystals, the phonon carried by electron is 
virtual when the energy of moving electron is lower than that of optical phonon which has 
a finite energy gap. Superficially for BEC system phonon is gapless, no matter how small 
the velocity of the impurity is, phonon is able to create. But it can also not be created for 
arbitrary small impurity velocity due to above stated BEC mechanics. So phonon number 
we calculate is also virtual. Corresponding behavior for phonon number is showed in Fig. 2 
and is very similar to the behavior of effective mass of impurity. This observation is consis- 
tent with our intuition about polaron: the more number of the phonon encompassing the 
impurity, the larger the force dragging the impurity and so the larger the impurity effective 
mass. 

Below we calculate energy related to an impurity with zero velocity. The whole energy 
of system is 

E g (0) = Eb + E(Q) + (j,n 
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ph 




(a) 



k B a , B 




Figure 2: The behavior of the phonon number Tipfi ELS ct function of boson-impurity scattering 
length aiB for different mass ratios /mi of boson to impurity (a); as a function of boson- 
impurity mass ratio tub /mi for different boson- impurity scattering length a/#. We have 
chosen clbb^b = 0.87. 



« ( 5 ng + £ B ) + (An + £(0)) (28) 

Where we have used ([5]) and make another approximation \xn ~ [in§ which is valid for dilute 
BEC. The first term on the right hand side in (|28h is the energy of pure bosonic system, 
so the second term gives impurity energy in BEC system. The calculation about energy 
must be careful. In Hamiltonian ([I]), we have used the contact interaction to represent 
the true potential for which its Fourier transformation should fall off at large momentum. 
This substitution leads to the ultraviolet divergence of the ground state energy so that 
a regularization must be forced. In fact, this divergence is not fundamental and can be 
regularized by way of absorbing more high order scattering contribution into scattering 
length [25] 



Ikclib^ 2m re dA 



A 



m red h 2 ^ k 2 



(29) 



After taking these measures, the impurity energy Ei after scaled by 27rfi, aBBn/ms is 

El = p{1 + m ) ^ + ^ [1 + m f ^sl. T1 m (30) 
mi k B a B B V3tt 3 ' 2 mi (k B a B B, ) 1 



where 



roc 

n (2) = / dyy 2 
Jo 



. [1 + (m B /mi)]yi ( y 2 + 2p)^ [{m B /mj)y + (y 2 + 2p)V2] 
The corresponding behavior of Ei is plotted in Fig. 3. 



(31) 
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Figure 3: The behavior of impurity energy Ej as a function of boson-impurity scattering 
length aiB for different mass ratio ms/mj of boson to impurity (a); as a function of boson- 
impurity mass ratio m^/m/ for different boson- impurity scattering length aj B - We have 
chosen a BB k B = 0.87. 



4 Discussion and Conclusion 

At last, we discuss the validity of our approximation (|12p in the frame of LLP theory. 
Following the method calculating n p h, (fTTj) equals to 

n = n + J2 4 + £(«fc + V D I /(*) I' " 2 E Wf{k)f{-k) (32) 

k k k 

After some predigestion and dividing both sides of (|32p by n 

1 - 8 „3/2/„ i, ^3/2 | 2 P n , m i^2 {k B a IB ) 2 n ( 3 ) 

P 2 m/ (1 + mB )2 (W 1 n (4) (33) 



2 v / 3vr 3 / 2 n m B m; (k B a B B) 1/2 P(l - rj) 



with 



= / dy- 



y 2 + P i 



y 2 + 2/3 [( y 2 + 2p )i/2 + (me/m/)?,] 2 - 4(m B /m/) 2 P 2 (l - r/) 2 

Jo y {y 2 + 2p)[^TYp + ^y} n| a 2/ + ^2+^_ 2 ap(i- ?) ) 

(34) 

Equation (|32H is the most basic equation from which we can study the effect of single 
impurity on condensate fraction. For dilute bosonic system, the condensate fraction p 
approaches unit. The density of a typical BEC system in experiments is 10 13 — 10 15 cm -3 . 
Although it is small in contrast to typical density of gas, liquid and solid, it remains a large 
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number. From ()33[) there is a factor n in the denominators of last two terms on the right 
hand side. If there are lots of impurities, the last two terms would also include a factor 
proportional to impurity density. So the last two terms are small and negligible in contrast 
to the first and second terms in the light of the existence of single impurity. In this way the 
approximation (|12|) is proven to be valid. 

In conclusion, by making Bogoliubov and forward-scattering approximations we have 
obtained an equivalent Hamiltonian which is similar to the Frohlich large polaron Hamilto- 
nian and calculate effective mass of impurity, the phonon number carried by the impurity 
and the impurity energy on a basis of LLP theory. In addition, we also point out the valid 
range of forward-scattering approximation that the velocity of the impurity must be smaller 
than Landau critical velocity of BEC system. 
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